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ABSTRACT 

We present spatial distribution of polycyclic aromatic hydrocarbons and ion- 
ized gas within the Galactic giant H II region NGC3603. Using the IRS in- 
strument on board the Spitzer Space Telescope, we study in particular the PAH 
emission features at ~5.7, 6.2, 7.7, 8.6, and 11.3 //m, and the [Ar n] 6.99 jum, 
[Ne il] 12.81 /an, [Ar ill] 8.99 /an, and [S iv] 10.51 jum forbidden emission lines. 
The observations probe both ionized regions and photodissociation regions. Sil- 
icate emission is detected close to the central cluster while silicate absorption is 
seen further away. We find no significant variation of the PAH ionization frac- 
tion across the whole region. The emission of very small grains lies closer to the 
central stellar cluster than emission of PAHs. The PAH/VSG ratio anticorre- 
lates with the hardness of the interstellar radiation field suggesting a destruction 
mechanism of the molecules within the ionized gas, as shown for low-metallicity 
galaxies by Madden et al. (2006). 

Subject headings: HII regions, ISM: individual: NGC 3603, ISM: dust, ISM: atoms, 
ISM: molecules, infrared: ISM, Telescopes: Spitzer 
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1. Introduction 

The investigation of star-formation feedback on the interstellar medium (ISM) in 
low-metallicity environments is of prime importance to understand and constrain the 
evolution of primordial systems. The most relevant diagnostics are drawn from the study of 
the various ISM components as a function of the physical conditions in starburst-dominated 
objects such as dwarf star-forming galaxies and giant H II regions. In this perspective, the 
mid-infared (MIR) domain, which gives access to both gas and dust components, provides 
one of the most powerful tools. 

Star-forming objects are characterized by strong MIR emission features attributed to 
polycyclic aromatic hydrocarbons (PAHs) in the photodissociation envelopes surrounding 
massive star clusters (Leger & Puget 1984; Allamandola et al. 1985, 1989; Puget & Leger 
1989; Tielens et al. 1999). The PAHs emitting in the 3-13/xm range contain from several 
tens up to several hundreds carbon atoms (see, e.g., Schutte et al. 1993). PAH molecules are 
mainly excited by far-UV radiation. The MIR emission features are due to the subsequent 
fluorescence of aromatic C— C (especially dominant at 6.2 and 7.7 fim) and peripheral C— H 
(3.3, 8.6, and 11.3 /im) fundamental vibrational and bending modes. 

The various PAH emission features are differently affected by the local physical 
conditions such as the hardness of the interstellar radiation field (ISRF), the dust 
temperature, or the PAH mixture (Hony et al. 2001; Peeters et al. 2002). It has been 
suggested that strong UV radiation is able to ionize the PAHs, while for weaker radiation, 
PAHs can be neutral or even negatively charged by accretion of a single electron (Bakes 
& Tielens 1994, 1998; Salama et al. 1996; Dartois & d'Hendecourt 1997). Laboratory 
experiments indicate that neutral PAHs show stronger C— H mode emission relatively to 
C— C modes while the inverse is true for ionized PAHs (Szczepanski & Vala 1993; Langhoff 
1996; Kim et al. 2001; Hudgins & Allamandolla 1999). As a result, the emission features 
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at 3.3 and 11.3 //m are thought to mainly originate from neutral PAHs while the emission 
features between 6 and 9 fim are due to ionized PAHs. Recently, it has been demonstrated 
that not only the UV radiation but also the metallicity and the dust extinction are able to 
significantly influence the PAH ionization fraction (Cox & Spaans 2006). 

Because of their chemical composition, PAH molecules are expected to be less abundant 
in low-metallicity environments. The PAH intensity in galaxies has been found to correlate 
with the ISM metallicity (see e.g., Madden et al. 2006; Engelbracht et al. 2005; Wu et 
al. 2006). This correlation can be due either to the low carbon abundance, to the harder 
ISRF from low-metallicity hot stellar atmospheres (see e.g., Schaller et al. 1992; Schaerer & 
Maeder 1992;Schaerer et al. 1993; Charbonnel et al. 1993), or to a combination of these two 
effects. Madden et al. (2006) showed that the radiation field itself has an impact on the 
PAH survival in various metallicity dwarf galaxies. 

The PAH spectrum is seen to vary not only from one object to another but also within 
a single object. Hence it is possible to investigate the variations of the PAH spectrum as a 
function of the physical conditions across a given region (see e.g., Joblin et al. 2000; Vermeij 
et al. 2002; Bregman & Temi 2005; Kassis et al. 2006). In a few objects, the PAH intensity 
has been found to decrease when the ISRF hardens (Verstraete et al. 1996; Madden et al. 
2006; Beirao et al. 2006), suggesting that PAH molecules are destroyed by high-energy 
photons. 

The Infrared Spectrograph (IRS; Houck et al. 2004) on board the Spitzer Space 
Telescope (Werner et al. 2004a) opened a new perspective in the extraction of small-scale 
regions within extended source. The wavelength range covered by the IRS gives the 
possibility to investigate the spatial distribution of the ionized gas, the molecular hydrogen 
gas, the PAHs, the silicate dust, and the very small grains (VSGs). One of the most 
interesting application enabled by the IRS is to understand the influence of the local 
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physical conditions, such as the hardness of the ISRF, on the molecular content. 

Galactic giant H II regions are ideal benchmarks for such a study. NGC 3603 is a giant 
H II region located ~7kpc from the Sun along the Galactic plane (Moffat 1983; Moffat et 
al. 1994; Drissen et al. 1995; Brandl et al. 1999). Oxygen abundance measurements range 
from 12 + log(0/H) = 8.39 to 8.52 (Melnick et al. 1989; Tapia et al. 2001; Garcia- Rojas 
et al. 2006) and imply a metallicity close to solar. More than 50 O and WR stars (Moffat 
et al. 1994) produce a Lyman continuum flux of 10 51 s~ x (Kennicutt 1984; Drissen et 
al. 1995), which is about 100 times the ionizing power of the Orion Trapezium cluster. 
Within the core of the cluster, the system HD 97950 contains several WR, 03, and many 
late O stars (Melnick et al. 1989; Hofmann et al. 1995). The massive stars influence the 
surrounding ISM morphology, notably by compressing the molecular clouds through stellar 
winds (Niirnberger & Stanke 2002). For this reason, the geometry consists in a complex 
arrangement of numerous, localized H II region — PDR transitions. Although its bolometric 
luminosity of Lboi ~ 10 7 L is only about one tenth of the luminosity of 30 Doradus, it looks 
remarkably similar to R136, the stellar core of 30 Doradus (Moffat et al. 1994). NGC 3603 
has often been referenced as the most massive, optically visible H II region in our Galaxy, 
due to its high luminosity and moderately low extinction of A v ~ 4.5 (Eisenhauer et al. 
1998), most of which is Galactic foreground extinction. 

In this paper, we investigate with the IRS the spatial variations of the MIR features 
across NGC 3603. After introducing the observations in $21 we describe the data analysis in 
§21 The MIR morphology is then investigated in £J4J The gas distribution and in particular 
the ISRF hardness is derived in §5j Distribution of dust and molecular features is discussed 
in §6j Finally, we question the PAH survival in §3, and apply our results to usual MIR 
diagnostic diagrams in §81 



- 6- 



2. Observations 

NGC3603 was observed as part of the IRS GTO program to study massive star 
formation in Local Group giant H II regions (PID63). Extensive imaging using the 
Infrared Array Camera (IRAC; Fazio et al. 2004) was also peformed (Brandl et al. 2007, in 
preparation), and a total of nine infrared bright sources (clusters, shock fronts, protostars, 
etc..) have been selected for follow-up spectroscopy with the high- and the low-resolution 
modules of the IRS (AORKEY 12080384). The complete spectral analysis will be discussed 
in Lebouteiller et al. (2007, in preparation). The long slits of the IRS low-resolution 
modules, however, cover a significant area of NGC 3603 and provide information on the 
spatial variations of spectral features across the region. 

We present an analysis of the positions #3 to #8 from the original follow-up 
observations (see the coordinates in Tableland the image of Figured]). Position $4 
is an offset observation and was used for sky substraction. Positions #2 and #9 (the 
latter corresponding to the IR luminous source IRS-9) gave corrupted and saturated data 
respectively. The position #6 is centered on the central stellar cluster while the other 
positions are centered on bright MIR knots located from ~1.6 to ~6pc away from the 
cluster (we use hereafter the projected distance from the central stellar cluster, which 
represents the smallest possible distance). 

The low-resolution spectra (A/AA ~ 65-130) from the short-low (SL) module cover 
the spectral region 5.2-14.5 /im. The SL module contains two slits, one for each spectral 
order. The SL1 slit corresponds to the first order (7.4-14.5 ^m) and has a size of 3.7" x 57" 
(<-» 0.13 x 1.93 pc 2 at a distance of 7kpc). The second order (5.2-7.7/xm) is observed 
through the SL2 slit (3.6" x 57" <-> 0.12 x 1.93 pc 2 ). The observations were done in staring 
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Fig. 1. — Slit positions are superimposed on an IRAC composite image. The colors are 
defined by the chl band ([3.17-3.93] /mi) for the red, chS ([5.02-6.44] /zm) for the green, and 
ch4 ([6.42-9.32] /zm) for the blue. For a given position, the northern slits observe only the 
second spectral order (SL2), the southern slits observe only the first spectral order (SL1), 
and the middle slit observations give access to both orders. The local PAH emission maxima 
(small circles) are plotted on the left panel. The size of the circles is proportional to the 
PAH peak intensity (sum of all bands). We also show the regions (arrows) showing an 
anticorrelation between PAH intensity and both [Ar ill] and [S iv] line intensities. 
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Table 1. Requested coordinates of each position. 



# 


a (J2000) 


5 (J2000) 


p A a 


l b 


ll h 14 m 21 s .80 


-61°24'20".0 


+24°.62 


T 


ll h 15 m 17 s .00 


-61°19'10".0 


+24°.42 


3 


ll h 14 m 49 s .06 


-61 o 17'09".l 


+24°.55 


4 


ll h 14 m 56 s .71 


-61°12'56".6 


+24°.52 


5 


ll h 14 m 52 s .40 


-61°15'46".3 


+24°.54 


6 


ll h 15 m 07 s .40 


-61°15'39".2 


+24°.48 


7 


ll h 15 m 08 s .03 


-61°16'40".2 


+24°.48 


8 


ll h 15 m 02 s .88 


-61°15'51".6 


+24°.50 


9 C 


ll h 15 m 03 s .30 


-61°21 / 25".0 


+24°.47 



a Position angle. 

b Position #1 was used for sky substraction. 
c Position #2 and #9 gave unusable data (see text). 
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modqj, consisting in two subsequent observations centered at the 1/3 (nod 1) and 2/3 (nod 
2) positions along the slits. When the source is being observed in a given slit (nominal 
observation), the other slit, corresponding to the other diffraction order, performs an offset 
observation. Since we are dealing with sources being more extended than the slit length, 
we take the opportunity given by the offset observations to extend our measures to larger 
spatial scales, ~ 2'. 



3. Data analysis 

3.1. Detector image reduction 

The data were processed by the Spitzer Science Center, using the S13.2 pipeline. We 
used the basic calibrated data (BCD) product. The background was calculated using the 
nominal and offset observations of position #1, and substracted from the other positions. 

Spectra were extracted from the two-dimensional flat-fielded detector image using 
scripts within the Cornell IRS Spectroscopy Modelling Analysis and Reduction Tool 
environment (SMART; Higdon et al. 2004). The two order spectra scale relatively well, and 
no corrections were needed. The full slit spectra of the various positions and the global 
spectrum of NGC 3603 (simple sum of the spectra of all positions) are shown in Figure[2] 
The final resolution is AA = 0.06 /im for SL2 and 0.12 /zm for SL1. 



3.2. Extraction of small-scale regions 

Spatial variations can be observed within the long slits (see Figure[3]), with regions 
emitting strong forbidden lines and no PAHs, and regions showing intense PAH features 



1 See the Spitzer/IRS observer's manual at http://ssc.spitzer.caltech.edu/documents/som/ 
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Fig. 2. — In the top panel, the full slit spectra of each position (nominal observations) are 
sorted as a function of the projected distance from the slit center to the central stellar cluster. 
In the bottom panel, we show the sum of these spectra. 




Fig. 3. — Spectral image of the nominal observation of position #3. The flux (light scale) 
is plotted as a function of the spatial position within the slit (y-axis) and of the wavelength 
(x-axis). The segment on the right side indicates the size of the extraction windows. 
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together with relatively weak lines. To investigate these spatial variations, we extracted 
columns, fixed to 4 pixels wide, within the long slit on the basis that the full width at half 
maximum of the point spread function (PSF) is ~2 pixels. MIR sources above or below 
the slit, or on each side of the extraction window within the slit, can contribute (because 
of the PSF) to the total flux in a given chosen extraction window. Hence the spatial 
column actually extracted corresponds to a region somewhat larger than 0.24x0.12 pc 2 — 
which would be the size of the extracted spatial region if the PSF width was null. For 
this reason, it is particularly difficult to estimate the related absolute uncertainty on the 
flux calibration. Given the fact that our study is focused on relative spatial variations, we 
consider this caveat as being of little importance as far as the interpretations are concerned. 
Comparing the fluxes of same spatial positions observed at different locations within the 
slit, we find relative errors ranging from ±7% (position #3, lowest integrated flux) to ±4% 
(position #7, largest flux). 

We extracted overlapping extraction windows along the slits, shifting by one pixel 
between each other (which make a total of 37 different columns, and as many data points, 
for a given nominal or offset observation). This is not an actual oversampling of the spatial 
PSF. Indeed, whatever the sampling is, one is limited by the PSF convolution. The spatial 
profile of spectral feature emission along the slit are smoothed by the PSF so that the choice 
of the number of extracting windows only results in different samplings of the — same — 
eventually smoothed spatial profile. Finally, it must be added that variations of features in 
the observed (PSF convolved) spatial profile imply in reality even larger variations. 

3.3. Measurements 

We observe several PAH features and forbidden emission lines superimposed on a 
spectral continuum which in our case is dominated by thermal emission of dust and broad 
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silicate emission/absorption. Several windows in positions #3, #4, #5, and #7 show 
silicate absorption around 10 /im while positions #6 and #8 show silicate emission in the 
same range (see Fig.H]). Note that the silicate features appear less prominent in the full 
slit spectra of Figure[2] due to the contribution of several distinct physical regions within 
the slit. Only the position #6 shows signs of a stellar emission continuum rising toward 
wavelengths shorter than ~8 /im (see Fig. [2]). 

PAH features are detected at ~5.7, ~6.2, ~7.7, ~8.6, ~11.3, and ~12.7//m. A weak 
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Fig. 4. — Examples of the spectral fit of a 4 pixels window in positions #3 (top) and #8 
(bottom). The various dust continua, PAH components, and forbidden emission lines can 
be identified. We also show the coverage of the IRAC/chS and ch4 bands. The chl and 
ch2 bands observe at wavelengths shorter than 5 /im. 
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feature at ~7.9/im can be observed for some positions as a shoulder of the 7.7/im PAH. 

The following forbidden emission lines are detected, [Ne n] (21.6 eV) at 12.81 /xm, 
[Ar II] (15.8 eV) at 6.99 /an, [Ar ill] (27.6 eV) at 8.99 /im, and [S iv] (34.8 eV) at 10.51 /im. 
The [Ne n] line is blended with the PAH feature at 12.7/xm. The [Ar n] line could possibly 
be blended with a weak PAH feature at 7.0 /im. This contamination is probably not larger 
than ?^20% of the total line flux (see the measurements of Sturm et al. 2000; Forster 
Schreiber et al. 2001). 

The intensity measurements were done using the spectral fitting algorithm PAHFIT 
(Smith et al. 2007). This procedure is especially suited for mid-IR low-resolution spectra 
of star-forming regions dominated by PAH emission. The modelled spectral continuum is 
computed from a combination of starlight and 8 preset (fixed temperature blackbodies) 
thermal dust continua, whose strength is allowed to vary. The unresolved forbidden 
emission lines are represented by Gaussian profiles, while the PAH features are decomposed 
in individual blended Drude profiles. The main components we consider in this study are 
the PAH bands at 5.70 /mi, 6.22 /im, 7.60+7.85 /im (which we will refer from now on as the 
7.7/im feature), 8.61 /im, and 11.23+11.33 /im (11.3 /im feature). The absorption by silicate 
dust, as opposed to its emission, is accounted for in the PAHFIT calculations. We decided 
to model the silicate emission in positions #6 and #8 by including two broad Gaussian 
components, centered at 9.5 and 10.5 /im. This is an arbitrary choice whose main purpose is 
to be able to measure the superimposed features with the smallest uncertainty possible. We 
integrated the PAHFIT algorithm in an automatic script able to handle the large number 
of spectra (>600) of our dataset. Examples of spectral fittings are presented in FigurelU 
The small-scale variations of the MIR spectral feature intensities across the slit positions 
are given in Table [2] and plotted in FigureEl 

Finally, in order to trace the presence of VSGs, which are expected to dominate the 
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Fig. 5. — Spatial variations of the spectral feature intensities. Integrated flux is plotted 
against the declination for each slit position. Data on the the left side of each plot correspond 
to the southern offset observation (first spectral order only, 5.2-7.7 /xm), and data on the right 
to the northern observation (second spectral order only, 7.4-14.5 /im). Legend is displayed 
in position #7. The bottom panels show the relative errors. The black rectangles along the 
declination axis indicate regions where PAH intensities spatially anticorrelate with [Ar ill] 
and [S iv] line intensities. The large vertical gray rectangles indicate zones with no coverage. 
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continuum emission longward of A ~ 10 /im (Desert et al. 1990), we measure the average 
flux between 13.6 and 14.0 /im, where there is no contamination by line emission and PAH 
features. 
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Table 2. Flux of the lines and the PAH features. 



# 


a, 5 (J2000) 




1 a 

d a 


5.7/xm 


6.2 /im 


7.7/xm 


8.6 /iin 


11.3 /iiii 




(degj 




(pc) 












3 


168.70210, - 


61.2883 


/"> f\f\f\ 

6.000 


0.2831 


1.9817 


4.2977 


1.0332 


1.5536 










(0.0516) 


(0.0739) 


(0 1 ^0^ 

IU. ±OVJO 1 


Cfl 0487^1 


fO 1 001 


Q 
O 


1 ftfi 709/1 /I 
lOo. U244, — 


ci noon 
Ol.ZOoU 


o.yoo 


0.3213 


2.1403 


A K7Q9 
4.0 1 oz 


l.zuio 


l.oU4z 










(0.0531) 


(0.0677) 


(n i qqi "\ 




(n nQ7rO 
^u.uy i u j 


• ) 


1 68 70644 - 


•61 2838 


5 466 


0.1318 


1.0722 


1 7679 


3879 


8580 










(0.0495) 


(0.0395) 


yU.L 1UU J 














(truncated) 








1 1 

# 


f" / Tl~»A AA \ 

a, 6 (J2000) 




1 a 

d a 


[Ar n] 


[Ne ll] 


r a TTT i 

[Ar III] 


[S iv] 






(deg) 




(pc) 












3 


168.70210,- 


61.2883 


6.000 


0.0350 


0.4403 


0.1572 


0.2307 












(0.0085) 


(0.0103) 


(0.0214) 


(0.0247) 




3 


168.70244, - 


61.2880 


5.958 


0.0371 


0.5342 


0.1861 


0.2610 












(0.0074) 


(0.0098) 


(0.0155) 


(0.0225) 




3 


168.70644, - 


•61.2838 


5.466 


0.0333 


0.6071 


0.2992 


0.5256 












(0.0071) 


(0.0063) 


(0.0150) 


(0.0220) 





(truncated) 



a Distance from the central stellar cluster (a = ll h 15 m 07 s .966, 5 = -61°15'30 .348). 
Note. — Fluxes are expressed in xlO 20 Wcm -2 . Numbers between brackets represent the 
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uncertainties. 
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4. Interpretation of the IRAC image 

The IRAC image of NGC3603 (Fig. [I]) reveals a complex MIR morphology that the 
spectroscopic results of the IRS in Figure[5] can help us to understand. 

We expect the IRAC/chl band to be dominated by stellar continuum emission and 
by the 3.3 /im PAH feature, which is seen to scale with the 11.3 /im feature in various 
objects (Hony et al. 2001). As an illustration, the full slit spectrum of the position #6, 
which is centered on the central stellar cluster, is strongly dominated by stellar emission at 
wavelengths shorter than 8/zm (Fig.[2|). Furthermore, Figured] shows that stars indeed emit 
mostly in the chl band. 

The regions where PAH intensity is the largest are bright in both chl and chS bands. 
This is because of the presence of PAH features (3.3 /im and 6.2 /zm) within these bands. 

The ch4 band is sensitive to the presence of the PAH features at 7.7 and 8.6 fim, and 
to the forbidden emission lines [Ar n] and [Ar ill]. By comparing the IRAC image with 
the results of the IRS, we notice that the regions showing a flux excess in the ch4 band as 
compared to the other bands are also those showing a relatively intense [Ar ill] line in the 
IRS spectra. 

5. Gas distribution 

The forbidden line intensities vary significantly across each observation. We assume on 
first approximation that the metallicity is uniform within the region and that it does not 
affect the spatial variations we observe. Preliminary results on abundance determinations in 
NGC3603 confirm the validity of this assumption (Lebouteiller et al. 2007, in preparation). 
The variations are neither due to depletion effects, at least for argon and neon which are 
not expected to be depleted onto dust grains. The depletion of sulfur is more uncertain. 
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Abundance determinations in H II regions and planetary nebula? suggest that sulfur 
depletion ranges from null to — 0.5dex with respect to the solar abundance (Pottasch & 
Bernard-Salas 2006; Henry et al. 2004; Marigo et al. 2003; Martin-Hernandez et al. 2002). 
On the other hand, the line intensity spatial variations and, to an even larger extent, the 
line ratios, are sensitive to the ionization structure. Since both high-ionization lines [Ar ill] 
and [S iv] follow the same trend, the depletion of sulfur should not dominate the spatial 
variations of [S iv] line intensity. 

Because the line ratios implying species with large differences in their ionization 
potentials can actually probe different regions along the sightline, we expect a priori the 
ratios [S lv]/[Ne il], [S iv]/[Ar ill], and [Ar m]/[Ne n] to be our most reliable tracers of 
the transversal spatial profile of the ISRF hardness. The usual MIR estimators used in 
the literature are [Ne lll]/[Ne n] and [S iv]/[S ill], which we cannot measure with our 
dataset. In order to compare the reliability of our estimators, we used the results of the 
stellar population and photoionization code developed by Guiles et al. (2004). They used 
Starburst 99 (Leitherer et al. 1999) and MAPPINGS Illq (Dopita et al. 2002; Sutherland 
et al. 2002) to derive the variation of the MIR photoionization lines as a function of the 
characteristics of the stellar population. We computed the line ratios in environments with 
metallicities between 1/3 and twice solar. The stellar population in the models was defined 
by an instantaneous burst, with a Salpeter initial mass function (Salpeter 1955), and an 
upper mass cut-off of 100 M Q . 

The models predict that the [S iv]/[S ill] ratio traces [Ne m]/[Ne n], as expected, 
and that the [S lv]/[Ne n], [S lv]/[Ar ill], and [Ar lll]/[Ne n] ratios correlate as tightly 
with [Ne m]/[Ne n]. They can be thus considered reliable tracers of the ISRF hardness in 
the environments considered. From now on, we will make use of [S Iv]/[Ne II], principally 
because this ratio involves the two most prominent forbidden emission lines in our spectra. 
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We plot in Figure[6] the distribution of [S iv] , [Ne n] , and their ratio as a function of 
the projected distance from the central cluster. It can be seen that I([S iv]) shows a sharp 
increase between and 1.5 pc, and then decreases progressively until distances larger than 
8 pc. The lack of emission at very small distances from the cluster is probably a consequence 
of the stellar winds that washed out the surrounding ISM to create a cavity. The image 
of Figured] shows indeed the total absence of interstellar MIR emission in the immediate 
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Fig. 6. — [S iv] (top), [Ne n] (middle) intensities, and the [S iv]/[Ne il] ratio (bottom) 
are plotted as a function of the projected distance from the central cluster. Intensities 
are expressed in xlO 20 Wm -2 . In order to minimize the projection effects, results in each 
distance bin of 0.5 pc were error-weighted and averaged. The vertical dotted lines refer to 
the distance marks plotted in Figure[U 
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vicinity of the central cluster. Similarly to [S iv], the [Ne il] line intensity is relatively 
weak at distances < 1 pc but then shows a fairly flat emission until ~6 pc. The shallow 
distribution of /([Ne il]) suggests that [Ne n] emission is more extended than that of [S iv] 
The [S lv]/[Ne il] ratio shows roughly the same behavior as the [S iv] line intensity, and 
implies that the ISRF hardness decreases by a factor ~20 toward the outer parts of the 
giant H II region. 



6. Dust and molecule distribution 
6.1. Silicate dust 

Amorphous silicate dust shows a broad spectral feature centered on 9.7/zm, originating 
from the stretching mode of the Si— O bond (see, e.g., Knacke & Thomson 1973). While it 
is mostly observed in absorption in astrophysical objects, silicate emission has also been 
detected in a few H II regions, including the Orion nebula, and is thought to be due to 
grains with size >0.01/im heated to ~100K (Cesarsky et al. 2000; Contursi et al. 2000). In 
NGC3603, we detect silicate dust in absorption in the spectra of positions #3, #4, #5, 
and #7. Silicate in emission is observed in the spectra of the other positions, #6 and #8. 
Examples of various silicate profiles across the region are presented in FigureO 

We measured the apparent strength of the 9.7/xm feature using the method of Spoon 
et al. (2007). This method requires to choose a local MIR continuum, and to calculate the 
ratio between the observed flux at 9.7/im and the continuum value. The silicate strength, 
as defined by the equation (1) of Spoon et al. (2007), is positive for silicate emission, 
and negative for silicate absorption. The spatial distribution of the silicate strength in 
NGC 3603 is presented in FigureO Silicate dust is seen in emission relatively close to the 
central cluster, while it is detected in absorption further away. 
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It is interesting to notice that the silicate emission is observed around the same 
distances to the cluster as the VSG emission ( §6.31) . The transition zone between silicate 
emission and silicate absorption is located between 2 and 3.5 pc away from the cluster. 
This corresponds to the region where the ISRF hardness decreases significantly (§3} and 
where the PAH emission begins to dominate the spectra. Finally, we observe that the most 
prominent silicate absorption features seem to correspond to bright PAH emission regions. 

6.2. PAHs 

For the following analysis, we consider hereafter that extinction has a negligible effect 
on the PAH spectrum from one position to another, or across a given observation (although 
it certainly has an effect on a single spectrum, see e.g., Peeters et al. 2002). 

n | i i i | i i n-| i i i i i i i n 




X (Hm) 



Fig. 7. — The 9.7/zm is detected in emission and in absorption. The two bottom spectra 
are from position #3, the two on the top from the position #8, and the middle one is from 
position #7. 
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The peak wavelength of the PAH features at 6.2, 7.7, and 8.6 /zm is expected to possibly 
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Fig. 8. — Silicate strength is plotted as a function of the distance from the central cluster. 
See FigureE] for the plot description. 
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Fig. 9. — (a) Intensity of the 11.3 /xm feature (neutral PAHs) is plotted against the sum of 
the 6.2, 7.7, and 8.6 /xm (cations) feature intensities, (b) The ratio 11.3/(6.2+7.7+8.6) is 
plotted against the distance to the stellar cluster. See FigureE] for the plot description. 
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shift, depending mostly on the molecule properties (size, symmetry, ...) as opposed to the 
11.3 /im profile which is found invariably peaking at the same wavelength (see, e.g., Hony et 
al. 2001). However, we find that the peak of all the PAH profiles does not vary by more than 
one resolution element across the observations of NGC3603 (0.06 /im for the 5.7 and 6.2 //m 
features, 0.12/im for the 8.6 and 11.3 /im features), the only exception being the PAH 
feature at 7.7 /im which shows a somewhat larger distribution, centered on 7.64±0.24/im. 
These results imply that the molecule properties do not change significantly across the 
region. In particular, the invariance of the 6.2 /im profile peak, centered on 6.23±0.06 /im, 
suggests the presence of complexed and/or substituted carriers (by introduction of an hetero 
atom in place of a carbon atom in the skeleton; Hudgins et al. 2005; Peeters et al. 2002). 
The spectra across NGC 3603 identify with the " class A" in the classification of Peeters et 
al. (2002), i.e., where the PAH features actually peak at r^6.22, «7.6, and ~8.6 /im. This 
class is populated by H II regions, reflection nebulae and most extragalactic sources. 

On first approximation, the PAH spectra across the region look identical, with the 
intensities of the various features scaling together. This suggests, together with the presence 
of both neutral (11.3 /im) and ionized (6.2, 7.7, 8.6 /im) features in the spectra, that the 
PAH ionization fraction is relatively constant. In order to investigate in more details the 
influence of ionization on the PAH spectra, we compare in Figured the intensity 7(11.3) 
with 7(6.2)+J(7.7)+J(8.6). There is a tight correlation, implying that the neutral/ionized 
mixture does not vary significantly across the region. More particularly, the ratio is 
essentially constant as a function of the distance to the cluster (Fig. [9b). We attribute the 
constant ionization fraction to the electron recombination rate in the ionized region which 
balances the PAH ionization rate. 

FigurefTUh shows that emission bands dominated by C— C modes and C— H modes 
correlate with each other. This correlation hold for the range of [S iv]/[Ne n] ratios probed 
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across the region ($2J). There is no correlation with the distance to the cluster (Fig.fTUb). 
These findings are consistent with the fact that the PAH ionization fraction is constant 
across NGC3603, since ionized PAHs should show enhanced C— C mode emission (see 
introduction). 



6.3. Very small grains 

While PAH emission originates in PDRs (Tielens 1993; Sellgren et al. 1990), VSG 
emission is seen to peak in ionized regions (Cesarsky et al. 1996; Verstraete et al. 1996). In 
order to check whether these results hold for NGC 3603, we estimated the VSG emission by 
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Fig. 10. — (a) The summed intensity of the 6.2 and 7.7/zm PAH features (C— C mode) is 
plotted against the sum of the 8.6 and 11.3 /xm feature (C— H mode) intensities, (b) The 
ratio (6.2+7.7)/(8.6+11.3) is plotted against the distance to the stellar cluster. See Figure[6] 
for the plot description. 
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measuring the continuum flux at 14 /xm. The VSG emission spans a wide range of values 
across the region (Fig. fTTj) . 

We find that the VSG continuum intensity scales tightly with [S iv] line intensity 
(Fig.lT2~l). Since [S iv] line intensity is seen to peak close to the central cluster (Fig.[6]), we 
conclude that VSGs also emit mostly in these regions and that they are the dominant dust 
component in the ionized region. 

7. PAH survival 

7.1. Comparison with VSG emission 

The difference between the spatial emission of PAH and VSG can be seen in Figure[T3| 
where we plot the PAH/VSG ratio as a function of the distance to the cluster. 

P I 1 ' ' I ' 1 ' I ' 1 ' I ' 1 r ~3 




Wavelength (flm) 

Fig. 11. — Example of various spectra across NGC3603, showing the diversity of the VSG 
continuum. 
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VSG and PAH emission do not coexist spatially, the VSGs emitting mostly in the 
ionized region (although they can be present elsewhere but without being excited). The 
photons exciting VSGs also illuminate PAHs, thus the PAH molecules may not survive the 
physical conditions required to heat VSGs (Madden et al. 2006). In order to investigate 
the relation between PAH emission and the hardness of the ISRF in more details, we use 
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Fig. 12. — /([S iv]), expressed in xlO 20 Wm 2 , is plotted against the VSG emission (in 

Jy). 
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Fig. 13. — The PAH/VSG ratio is plotted against the distance to the central cluster. See 
Figure[6]for the plot description. We use the PAH features at 8.6 and 11.3 /xm only since we 
plot here the results of the nominal and SL1 offset observations. 
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a similar approach to that of Madden et al. (2006) by comparing the PAH emission with 
VSG emission. We show in Figured!] the variation of the PAH/VSG intensity ratio as a 
function of [S lv]/[Ne ii]. It can be seen that PAH emission becomes globally weaker when 
the ISRF becomes harder. 

The PAH features dominated by C-C or C-H emission modes do not show significant 
differences across the region ( §6.2)) . Thus, dehydrogenation (rupture of the C— H bonds by 
absorption of UV photons) is unlikely to be responsible for the lack of PAHs relative to 
dust in regions where the ISRF is the hardest. On the other hand, PAH molecules could be 
destroyed by high-energy photons. 



We find that [Ar n] and [Ne n] emission peaks always correlates with PAH emission 
peaks. We find no regions where PAH emission anticorrelates with both I ([At ii]) and 
/([Ne ii]). In contrast, some regions show an anticorrelation between PAH intensity and 
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Fig. 14. — The PAH intensity (sum of the 8.6, and 11.3 /zm features) is plotted against the 
ISRF hardness probed by the [S iv]/[Ne ii] ratio. 
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J([Ar ill]) and 7([S iv]) (see their identification in Fig.[5]and their corresponding locations 
in Fig.[T]). Within these regions, the location of the ion peak emission always follows a clear 
structure, being correlated with the ionization potential, and implying a sharp variation of 
the hardness of the ISRF. 

We identify these transitions as being interfaces between ionized region and PDRs. The 
interfaces toward the ionized region regime are all located in the direction the central cluster 
(Fig. [I]). The size of the interfaces (between the maximal and minimal PAH intensities) is 
respectively 0.48, 0.72, 0.57, 0.66, and 0.40 pc (note that this should not be confused with 
the size of the PDR). 

7.3. Photodestruction 

The energy deposited via photon heating is potentially large enough to dissociate 
bonds within PAHs. This has been the explanation to the PAH variations in the Orion bar 
and in the Red Rectangle nebula (Geballe et al. 1989). The Ar II ion exists for energies 
larger than 15.8 eV and it is the dominant ionization stage in regions where PAH intensity 
is maximal. This means that the far-UV radiation responsible for the dominant presence of 
Ar II (from 15.8 to 27.6 eV which is the ionization potential of Ar il) is not able to destroy 
efficiently the PAHs, at least for the intensity of the ISRF at this energy. Although the 
energy required to break C— H or C— C bonds is 4.8 eV and ~8.5eV respectively (Omont 
1986), photons with higher energy may be needed to dissociate bonds in a large molecule. 
Using the models of Omont (1986), we find that bonds can be broken by < 27.6 eV photons 
only in a molecule smaller than ~25-50 C-atoms, depending on the bond type. Note that 
the threshold size varies from one model to another, and the values we derive should only 
give a first order approximation. Hence, one possibility is that PAH molecules are large 
enough to prevent dissociation to occur. Another possibility requires that PAH molecules 
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are small, and that the energy density is relatively weak. 

Can we infer a maximal size for the PAHs in NGC 3603? The UV radiation required 
for the dominant presence of Ar ill (27.6-40.7eV) and S IV (34.8-47.3 eV) ions could be 
responsible for PAH molecule destruction. Single photons with energies 40.7 eV (47.3 eV) 
are able to break bonds in molecules smaller than ~40-75 (~50-85) C-atoms, depending 
on the bond type. Furthermore, it also implies that the energy density is large enough. 
Although even larger PAHs may be present, it is difficult to set an upper limit on their size. 
However, since PAH emission is almost zero in regions where Ar ill and S IV are dominant, 
it is likely that most PAHs have sizes smaller than ~85 C-atoms, at least in regions where 
the ISRF is the hardest. 

8. MIR diagnostic diagram 

The small-scale variations of the spectral feature intensities across NGC 3603 show 
the transition between regions dominated by PAH emission and regions dominated by 
high-ionization lines. Our results should fit in MIR diagnostic diagrams used to distinguish 
between active galactic nuclei (AGN), PDR, and ionized region behaviors. 

In Figure[15l we present a diagram, based on Laurent et al. (2000) and Peeters et al. 
(2004). The templates we use are the Galactic reflection nebula NGC7023 (Werner et al. 
2004b; PDR excited by optical photons), the H II region M17 - position #2 (Peeters et al. 
2004), and the quasar 3C273. Accordingly to the approach of Laurent et al. (2000) and 
Peeters et al. (2004), we plot the continuum ratio between ~i3.8/im (average flux between 
13.6 and 14.0/zm) and ~6.0 fim (average between 5.9 and 6.1 /im) against the ratio of the 
6.2 /iin PAH intensity over the continuum level at ~6.0 /mi. 

Our data points form a relatively narrow stripe, probing observationally the ionized 
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Fig. 15. — MIR diagnostic diagram. The dashed and dotted curves indicate 50-50% and 25- 
75% contributions respectively. Small filled circles correspond to individual column spectra. 
The stars refer to the columns within position ^6 which are centered on the central stellar 
cluster. Large filled circles represent the measure using full slit extractions. The diamond 
stands for the measure of the global summed spectrum of NGC3603. 
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region - PDR transition (Fig. [T5|) . The results of positions #7 and #8 lie around the ionized 
region template, while the results of positions #3, #4, and #5 range from 100% ionized 
region-contribution to 75% PDR-contribution. There is no region within our observed field 
in NGC3603 similar to an isolated "exposed" PDRs such as NGC7023. This is likely due to 
a geometry effect because dust emission at 14 /im could lie in the foreground or background, 
resulting in an overestimation of the 14.0/6.0 continuum ratio. 

The only points showing more than 50% "AGN-like" behavior are those from position 
#6 that are probing the central stellar cluster. This is because for these positions, the 
continuum measured at 6 /im is dominated by stellar emission. Thus, if a given starburst 
galaxy spectrum is dominated by stellar cluster emission, it could a priori be confused with 
an AGN when using MIR diagnostics. However, the global spectrum of NGC 3603, obtained 
using the sum of the full slit spectra of all positions (diamond in the Fig. [T5"j) . does not show 
any sign of stellar contamination. In this case, the regions characterized by dust emission 
dominate the global spectrum. Other than the small-scale results of position #6, none of 
the points probing purely ISM material show signs of AGN regime. Using the observations 
in NGC 3603, we successfully test this diagram as a diagnostic tool to distinguish regimes 
in a single objects, as long as the stellar emission is not significant. 

We also show in FigurefTSl the results using full slit extraction of each position (large 
squares), that were obtained by integrating the MIR emission of all the regions within the 
slit. The corresponding data points of each position lie in the middle of the small-scale 
results, implying that the full slit extraction does not give any systematic error on the 
diagnostic. 



-33- 



9. Conclusions 

We have investigated the spatial variations of the MIR spectral feature intensities 
(both ionic and solid state) in the Galactic H II region complex NGC3603, and obtained 
the following results: 

1. On first approximation, the various PAH emission features have identical spatial 
distribution and the PAH ionization fraction is constant. 

2. The ISRF hardness, as probed by the [S iv]/[Ne n] ratio, decreases as a function of 
the distance from the central stellar cluster. 

3. Silicate is detected in emission close to the cluster while it is detected in absorption 
further away. Local absorption maxima of the 9. 7 /mi feature seem to identify with 
bright PAH emission knots. 

5. PAH emission lies at larger distances than VSG emission and becomes weaker 
when the ISRF becomes harder. It has been shown for low-metallicity galaxies by 
Madden et al. (2006). It seems that PAH molecules are not able to survive the same 
physical conditions as VSGs. The simplest explanation is the photodestruction of the 
molecules. 

6. Small-scale results within NGC3603 allowed us to probe observationnally the 
transition between ionized region and PDR regimes in the MIR diagnostic diagram 
of Laurent et al. (2000). In this diagram, the measurements using individual column 
spectra form a relatively narrow stripe. 

We thank E. Peeters for providing data for the diagnostic diagram (Fig. [T5|) . S. Madden, 
G. Sloan, and H. Spoon for their useful comments. This work is based on observations 
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California Institute of Technology, under NASA contract 1047. Support for this work was 
provided by NASA through contract 1257184 issued by JPL/Caltech. 
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